TORIC GEOMETRY OF CUTS AND SPLITS 



BERND STURMFELS AND SETH SULLIVANT 

Abstract. Associated to any graph is a toric ideal whose generators record relations 
among the cuts of the graph. We study these ideals and the geometry of the correspond- 
ing toric varieties. Our theorems and conjectures relate the combinatorial structure of 
the graph and the corresponding cut polytope to algebraic properties of the ideal. Cut 
ideals generalize toric ideals arising in phylogenctics and the study of contingency tables. 



1. Introduction 

With any finite graph G = {V, E) we associate a projective toric variety Xq over a 
field IK as follows. The coordinates qA\B of the ambient projective space are indexed 
by unordered partitions A\B oi the vertex set V . The dense torus has two coordinates 
(sjj, tij) for each edge j} G E. The polynomial rings in these two sets of unknowns are 

K[g] := K[qA\B\ AU B = V,Ar}B = 
K[s,t] := K.[sij,Uj\{i,3}^E]. 

Each partition A\B oi the vertex set V defines a subset G\xt{A\B) of the edge set E. 
Namely, C\xt{A\B) is the set of edges {«, j} such that iEA, jEBoTiEB,jEA. The 
variety we wish to study is specified by the following homomorphism of polynomial rings: 

(1.1) : ^ qA\B ^ n ■ n 

{i,j}eCnt{A\B) {i,j}£E\Cnt(A\B) 

One may wish to think of s.. and t.. as abbreviations for "separated" and "together". The 
kernel of 0g is a homogeneous toric ideal Ig which we call the cut ideal of the graph G. 
We are interested in the projective toric variety Xq which is defined by the cut ideal Iq- 

Example 1.1. Let G = be the complete graph on four nodes, so V = {1, 2, 3, 4} and 
E = {12, 13, 14, 23, 24, 34}. The ring map (j)K^ is specified by 

'?|1234 ^ ^12^13^14^23^24^34! 9l|234 ^ ■512^13514^23^24^34! 

Q'l2|34 ^ ^12Sl3Si4S23S24^34! Q'2|134 ^ ■5i2tl3tl4S23'524^34! 

9l3|24 ^ ■Si2tl3Si4S23^24S34! 5'3|124 ^ ^12'5l3tl4S23^24S34! 

Q'l4|23 ^ Si2Si3ti4t23'524'534! Q'4|123 ^ ^12^13'5l4i23'524'534- 

The cut ideal for the complete graph on four nodes is the principal ideal 

= ( Q'|1234 '?12|34 Q'l3|24 Q'l4|23 — '?1|234 Q'2|134 Q'3|124 Q'l23|4 ) • 

Thus the toric variety Xk^ defined by Ik^ is a quartic hypersurface in P"^. □ 
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Example 1.2. Let G = C4 be the 4-cycle with edges E = {12, 23, 34, 14}. The ring map 
is derived from 0^-4 in Example 11.11 by setting S13 = tis = S24 = ^24 = 1, and we find 

= ( ?|1234 Q'l3|24 ~ 911234 Q'l24|3; Q'|1234 Q'l3|24 ~ 9l23|4 9l34|2; Q'|1234 Q'l3|24 ~ Q'l2|34 Q'l4|23 ) • 

Thus the toric variety is a complete intersection of three quadrics in P^. □ 

We usually take the vertex set V of our graph G to be [n] := {1,2,..., n}, so that K[q] 
is a polynomial ring in 2"~^ unknowns, and K[s, t] is a polynomial ring in 2\E\ < n(n — 1) 
unknowns. Each edge {i,j} G E corresponds to a projective line with homogeneous 
coordinates {sij : tij), and the ring map (pc represents a rational map from the product 
of projective lines (P^)'^' into the high-dimensional projective space P2"^ -1. The image 
of this map is our toric variety Xq, which has dimension \E\ < n{n — l)/2 in P^" 

The algebraic properties of its ideal Iq and the geometry of Xq are determined by the 
cut polytope Cut'~'(G), which is the convex hull in M'^' of the cut semimetrics Sa\b- Here 
A\B runs over all unordered partitions of V, and Sa\b ^ {0, 1}'^' is defined as follows: 

SA\B{{i,j}) = 1 if \Ar]{i,j}\ = 1 and SA\B{{i,j}) = otherwise. 

Indeed, the convex hull of the exponent vectors in (pc is affinely isomorphic to Cut'~'(G). 
In Example 11.11 and II. 2[ we find that Cut'~'(i^4) is the cyclic 6-polytope with 8 vertices, 
and Cut^(C4) is the 4-dimensional crosspolytope (which is the dual to the 4-cube). 

The cut polytope Cut^(G) is well-studied in combinatorial optimization, and is a central 
player in the book Geometry of Guts and Metrics by Deza and Laurent [7]. The title of 
this paper is a reference to their book, and refiects our desire to import this body of 
work into commutative algebra and algebraic statistics. In particular, we explore the 
extent to which the known polyhedral structure of Cut^(G') can be used to determine 
algebraic results about the cut ideals Iq- For instance, the known fact that Cut^(G') 
is full-dimensional implies that dimXc = \E\. A more significant example of such an 
algebraic result is derived from recent work of the second author [20] : 

Theorem 1.3. The cut ideal Iq has a squarefree reverse lexicographic initial ideal if and 
only if the graph G is free of minors and every induced cycle in G has length three or 
four. In this case, every reverse lexicographic initial ideal of Iq is squarefree. 

Proof. The initial ideal of a toric ideal is squarefree if and only if the corresponding regular 
triangulation of the associated polytope is unimodular [TSl §8] . Since the symmetry group 
of Cut^(G) is transitive on its vertices, the cut polytope Cut^(G) has a unimodular revlex 
(pulling) triangulation if and only if every revlex triangulation is unimodular [20l Cor. 2.5]. 
A polytope all of whose revlex triangulations are unimodular is called compressed. Now 
simply apply the classification of compressed cut polytopes given in [201 Thm. 3.2]. □ 

The outline of this paper is as follows. In Section 2, we describe how generating sets 
(Markov bases) and Grobner bases of the cut ideal Iq can be computed when the graph 
G admits a certain clique-sum decomposition. The key tool here is the toric fiber product 
which was introduced in [5T]. In Section [31 we summarize the results of our computational 
experiments, and we outline some conjectures which were suggested by our computations. 

In the last two sections we present applications to algebraic statistics. In Section [4] 
we relate cut ideals to the binary graph models of [6] and to Markov random fields. In 
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Section [5] we relate cut ideals to phylogenetic models on split systems, due to Bryant [2] . 
These generalize the binary Jukes-Cantor models which were studied in [1] and |19j . 

2. Clique Sums and Toric Fiber Products 

Our goal in this section is to relate the graph-theoretic operation of taking clique sums 
to the ideal-theoretic operation of taking the toric fiber product, as explained in [2l] . 
This operation will serve as a tool for reducing the computation of the cut ideals Iq to 
cut ideals of smaller graphs (and that, hence, involve fewer indeterminates) . 

Let Gi = (yi,Ei) and G2 = (V2,-E2) be graphs such that Vi fl V2 is a clique of both 
graphs. The new graph G = G'i#G2 with vertex set V = V1UV2 and edge set E = E1UE2 
is called the clique sum of Gi and G2 along fl V2. If the cardinality of Vi fl is k+1, 
this operation is also called a k-sum of the graphs. We suppose throughout that k <2. 

We now explain how binomials in the cut ideal Iq can be constructed from binomials 
in the smaller ideals Id and Jca- Consider an arbitrary binomial of degree d in the first 
smaller cut ideal Jci, say 

d d 
i=l i=l 

Since Vi fl V2 is a clique in Gi of cardinality < 3, we can permute the unknowns and 
partitions so that fl Vi fl V2 = Cj fl Vi fl V2 for all i. This is a consequence of the fact 
that Ir^j^^ is the zero ideal for k < 2. For any ordered list EF of d partitions of V2\Vi, 

EF = [Ei\Fi, E2\F2, . . . , Ed\Fd), 
we define a new binomial which is easily seen to be in the cut ideal Iq of the big graph: 

d d 
^EF _ ]^q'A,U£;,|B,UF, - J^Q'CUE.IDiUi^r 

i=l i=l 

This construction works verbatim if we switch the components Gi and G2, so that, for 
any binomial f in and any ordered list EF of deg(f) partitions of Vi\V2, we get a 
binomial f^^ in Iq- Moreover, if F is any set of binomials in or in Iq^ then we define 

(2.1) Lift(F) := { f^^ I f e F , = {Ei\F,}t\' } 

as the union of all binomials of the form f^^ described above. 

We also define an additional set Quad(Gi,G2) of quadratic binomials in Iq as follows. 
Let A\B be any unordered partition of Vi fl V2, let Gi\Di and Ei\Fi be any ordered 
partitions of Vi\V2, and let C2I-D2 and -E'2|-^2 be any ordered partitions of V2\Vi. Then 

(2.2) (lAUCiUC2\BUDiUD2 ■ QAUEiUE2\BUFiUF2 ~ (lAUEiUC2\BUFiUD2 ' (lAUCiUE2\BUDiUF2 

is in Quad(Gi, G2), and these are all the binomials in Quad((j'i, G2). For each fixed A\B, 
we can express the quadrics (12.21) as the 2 x 2-minors of a certain matrix {q,\,) of format 
2l^2\v^i| X 2l^i\^2|^ The following theorem will be our main result in Section 2. 

Theorem 2.1. Let G = G'i#G'2 0, 0, 1, or 2-sum of Gi and G2 and suppose that Fi 
and F2 are binomial generating sets for the smaller cut ideals Iqi and Ig2- Then 

M = Lift(Fi) U Lift(F2) U Quad(G'i,G'2) 
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is a generating set for the big cut ideal Iq- Furthermore, if Fi and F2 are Grobner bases 
for and then there exists a term order such that M zs a Grobner basis for Iq- 

Remark 2.2. If the intersection graph Gi r\G2 is not a chque of cardinahty < 3, then it 
is generally not possible to lift every binomial in and 7^2 to the cut ideal Iq- 

Before presenting the proof of Theorem 12 .11 we discuss several examples and corollaries. 

Example 2.3. If G = G'i#G'2 is a zero sum, then its cut ideal Iq is the usual Segre 
product of and Jgj. Indeed, in this case the singleton Vi fl V2 has only one ordered 
partition and Quad(Gi, G2) is the ideal of 2 x 2-minors of the corresponding matrix [q,],)- 
For instance, if Gi is the graph with one edge {1,2} and G2 is the graph with one edge 
{2, 3}, so 1^1 n V2 = {2}, then Iq = (Quad(Gi, G2)) is generated by the determinant of 

Now suppose that G is any tree with n leaves. Iterating the zero sum construction from 
77, = 3 to n > 3, we see that Xq is the Segre embedding of (P^)""-"^ into P^" □ 

Further generalization of Example 12.31 leads to the following result. 

Corollary 2.4. The toric variety Xq is smooth if and only if G is free of G4 minors. 

Proof. We first prove the if-direction. If G is free of C4 minors, so all its simple cycles 
have length three, then G can be built from K2 and by taking repeated 0-sums. Both 
the ideals and Ik^ are zero and live in polynomial rings with two and four unknowns 
respectively. Thus is P^ and X^^ is P^. The 0-sum construction amounts to taking 
Segre products, hence 

Xg = P^ X P^ X ■ ■ • X P^ X P^ X P^ X ■ ■ ■ X Pi 

This Segre variety is smooth. The only-if direction says that any smooth Xg has this 
special form. To prove this, suppose that G has C4 as a minor. Then either G has an 
induced cycle of length n > 4, or G has as an induced subgraph the complete graph 
K4 or the graph which is obtained from K4 by removing one edge. Let H denote this 
induced subgraph. Using a forward reference to Lemma 13. 2^ we note that Cut^(iJ) is 
a face of Cut'~'(G). Therefore, it suffices to prove that Xh is not smooth. We saw 
in the Introduction that Cut'~'(ii'4) and Cut^(C4) are not simple. Using the familiar 
characterization of toric singularities [9l §2.1], this implies that the corresponding toric 
varieties Xh are not smooth. The same can be checked for cycles of length n > 5. 

In the remaining case, H = i^'4\{14} is the 1-sum of the triangle on {1,2,3} and the 
triangle on {2,3,4}. Its variety X^ is the complete intersection of two quadrics in P^: 

Ih = ( det ( ^'^'A , det ( ^'^''^ ^''^A ) . 

\ \Q'4|123 Q'l4|23/ \Q'l3|24 93|124/ / 

The singular locus of Xh consists of the two 3-planes in P^ where these matrices are zero. 
The cut polytope Cut^{H) is the free join of two squares, a non-simple 5-polytope. □ 

The following example naturally generalizes the graph H = K4\{14} we just discussed. 
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Example 2.5. Let G = K5\{15} be the graph on five vertices obtained from the complete 
graph by deleting an edge. Thus G is the 2-sum of the complete graph Gi on Vi = 
{1,2,3,4} and the complete graph G2 on V2 = {2,3,4,5}. Since J/^^ is generated by a 
quartic, we deduce that Iq is generated by quadrics and quartics. There are four quadrics: 

Quad(Gi,G2) = {q'i5|234 9112345 ~ Q'112345 Q'511234 , Q'341125 9211345 ~ 9121345 9251134 

924| 135 93| 1245 ~ 9l3|245 935|124 , 923|145 94|1235 ~ 9l4|235 945|123 } • 

The ideals Iq^ and are each generated by a single quartic, as in Example II. ![ and 
Fi and F2 are the singletons consisting of these quartics. Now, the set V2\Vi = {5} has 
two ordered partitions, namely 5| and |5, so there are 2^ = 16 ordered lists of ordered 
partitions E\F. Each defines a quartic in 1^, so Lift(Fi) consists of 16 quartics, such as 

fl = 9|12345 934|125 924|135 923|145 ~ 9l|2345 925|134 935|124 945|1235 
h = 95|1234 9l2|345 9l3|245 9l4|235 ^ 9l5|234 92|1345 93|1245 94|1235- 

Likewise, Lift(F2) consists of 16 quartics, and these include 

f3 = 9l|2345 925|134 935|124 945|123 ~ 9l5|234 92|1345 93|1245 94|1235) 
f4 = 9|12345 934|125 924|135 923|145 ^ 95|1234 9l2|345 9l3|245 9l4|235- 

We conclude that the set M in Theorem 12. II consists of 36 binomials, and these binomials 
generate Iq- However, they are not a minimal generating set, because of the relation 

fl - f2 + fa - f4 = 0. 

The set of 35 binomials obtained by removing any of the fj is a minimal generating set for 
Ig- We also find that the minimal free resolution of Iq has the following Betti diagram: 

total: 1 35 134 200 134 35 1 

0: 1 

1: . 4 

2: . . 6 . ... 

3: . 31 128 200 128 31 . (Macaulay 2 output). 

4: . . . . 6 . . 

5: 4 . 

6: 1 

Thus the toric 9- fold Xq C P^^ is arithmetically Gorenstein. The degree of Xg is 80. □ 

We shall derive Theorem 12.11 from the results in [5T]. Specifically, we shall identify the 
cut ideal of G = Gi4^G2 as a toric fiber product. We begin by reviewing the set-up of 
[2Tj . Let r > be a positive integer and s, t G N*" be two vectors of positive integers. Let 

K[x] =K[xi\te [r], J e [s,] ] and K[y] = K[yi \ t e [r],k G [U] ] 

be two polynomial rings which have a compatible d-dimensional multigrading 

deg(x}) = degivl) = e (for ^ = 1, 2, . . . , r). 

We abbreviate the collection of degree vectors by ^ = {a^, a^, . . . , a^} C Z*^. 
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If / and J are homogeneous ideals of K[x] and ]K[?/] respectively, then the quotient rings 
R = K.[x]/I and S = K.[y]/J are also multigraded by A. Consider the polynomial ring 

K[z] = K[4 \ie[r],j e[s.^,ke[t,]], 

and consider the K-algebra homomorphism 

: K[z] ^R^kS, ^ x\ ® yl 
The kernel of is called the toric fiber product of / and J and is denoted 

/ J = ker(0/,j). 
The following statement combines Theorem 2.8 and Corollary 2.10 in |21j . 

Theorem 2.6. Suppose that the set A of degree vectors is linearly independent. Let Fi be 
a homogeneous generating set for I and F2 be a homogeneous generating set for J. Then 

M = Lift(Fi) U Lift(F2) U Quad,^ 

is a homogeneous generating set for I J . Furthermore, if Fi anc? F2 are Grobner bases 
for I and J, then there exists a term order such that ISA is a Grobner basis for I J. 

Here Quad_4 is the collection of quadrics Zj/^z}^ — Zj^^zl/^ which generates (0) x_4 (0). 
The sets Lift(Fj) have a nice description in terms of tableaux which is given in [211 §2]- 

Proof of Theorem \2.1[ Suppose G = Gi#G2 with vertex set V = V1UV2 and edge set 
E = E1UE2 where Vi fl V2 is a clique of size /c + 1 in both graphs. We set d = C^^^) + 1 and 
r = 2^~^, and we define A as the vector configuration corresponding to the vertices of the 
cut polytope Cut'~'(-ft'fc+i) of the clique. The ^-grading on K[q] is defined by restricting 
the product in (11.11) to those edges {i,j} which lie in Eir\E2. In other words, the degree of 
qA\B is the vertex of Cut°(i^fc+i) which is indexed by the partition A fl Vi fl V2 |-B fl Vi fl V2. 

The configuration A of degree vectors is linearly independent if and only if the cut 
polytope Cut^(i^fc+i) is a simplex if and only if k < 2. Theorem 12.61 requires the set A 
to be linear independent. This explains the crucial hypothesis /c < 2 in Theorem 12.11 

All three cut ideals Iq, Igi and 1^2 are homogeneous with respect to the indicated 
grading. We will show that Iq is the toric fiber product of 1^ and Jgj, in symbols, 

(2.3) Ig = Iq, x^Ig,. 

Let and A2I-B2 be partitions of Vi and V2 such that deg(gAi|Bi) = deg(gyi2|B2)- Since 

VLnV2 is connected, this implies (possibly after relabehng) that ylinVinV2 = A2nVinl^. 
This means that A\B with A = AiU A2 and B = BiU B2 is a. partition of V, and we have 

(2.4) 0Gi(q'Ai|Bi) ■ 0G2(?yi2|B2) = (pciqAlB) ■ <pGinG2iQAinVinV2\BinVinV2)- 

This is an identity of monomials in the polynomial ring K[s,t] associated with the big 
graph G, and it is verified by plugging in the definition of the monomial map 0, in (11.11) . 
The ring map which defines the toric fiber product Ig^ xaIg2 can be written as follows: 

(I^Igi,Ig2 ■ ^fe] ^[s^t]^ UlB ^ 0Gi(gAi|i?i) ■ 0G2(9A2|i?2)• 
Since (12. 4p holds and since 0GinG2(?AinViny2|BinViny2) divides (f)GiqA\B), the unknowns Sij 
or tij with {i,j} G Ei\~\E2 can appear in 0Gi (^AilsJ ■0G2(9a2|-B2) only with exponent 2. If 
we replace these unknowns Sjj, tij by their square roots in the monomial map (pj^ then 



TORIC GEOMETRY OF CUTS AND SPLITS 



7 



the kernel remains unchanged, and we get the monomial map (pc '■ K[s,t]. We 

conclude that ker(0G) = ker(0/g^^/g^), which is our claim (12.31) . Since the configuration 
A is linearly independent, we have thus derived Theorem 12. II from Theorem 12. 6[ □ 

The proof of Theorem 12.61 given in [21] reveals the possible choices of term orders which 
create a Grobner basis for Iq from given Grobner bases Fi of Iq^ and F2 of Jgj. First of 
all, the passage from a binomial f in Fj to the corresponding binomials f^^ in Lift(Fj) is 
compatible with the choice of leading terms, that is, we declare the leading term of to 
be the one coming from the leading term of f . In this manner we specify a family of partial 
term orders on K[g]. We then choose any tie-breaking term order on K[g] which makes 
the set Quad(Gi,G2) into a Grobner basis. Since these quadrics are the 2 x 2-minors of 
matrices {q,\m) whose entries are disjoint sets of unknowns there are many such choices 
of term orders. Now, the term order on K[g] which is gotten by refining the partial term 
order by the tie-breaker has the desired property that M is a Grobner basis for Iq- 

3. Computations and Conjectures 

Upon encountering a new family of ideals, our first instinct is to use computer algebra 
to gain a better "feel" for the way the structure of the ideals depends on the parameters 
defining the ideals. The parameter for the cut ideal Iq is the graph G, and we are 
interested in how the combinatorial structure of G determines the algebraic structure of 
Ig- To this end, we undertook an exploration of the cut ideals by computing generating 
sets, Grobner bases, free resolutions, and normalizations, using the programs 4ti2 [TT] . 
CoCoA [5], Macaulay 2 [TU], and Normaliz [T]. In this section, we summarize the results 
of our computations, and we offer a number of conjectures that arise from looking at the 
resulting data. 

3.1. Computations. The results are summarized in Table 1 below. The first column 
lists the graphs which we analyzed. These were all graphs on six or fewer vertices that 
are not clique-sum decomposable with a clique of size < 3. The notation of the form Gk 
comes from the Atlas of Graphs [TB]. However, if a graph has a more standard shorthand, 
we preferred to use the more easily identifiable abbreviations. The notations we used are: 

• Ki Complete graph, 

• Ki^^^^^^i^ Complete m-partite graph, 

• Ci Cycle of length /, 

• G Suspension of G over a point, 

• G X H Cartesian product graph. 

The columns in the table list the following features of the cut ideal Iq'- 
2-6 Number of minimal generators of Ig in degrees 2, 4, 6, 8, and 10. 

7 fi^Ic) = Largest degree of a minimal generator of Ig- 

8 Codimension (height) of Ig- 

9 Projective dimension of Ig- 

10 Degree (multiplicity) of Ig- 

11 Whether the semigroup algebra K[g]//G is normal. 

12 Whether the semigroup algebra K[q]/lG is Cohen-Macaulay. 

13 Whether the semigroup algebra K[g]/JGr is Gorenstein. 
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Table 1. Algebraic properties of cut ideals Iq for graphs G with up to six vertices. 



Blank spots in the table are entries that we were unable to compute. 

If G is a small clique-sum decomposable graph then we can break it into pieces that 
are listed in Table 1. This tells us the degrees of the minimal generators of cut ideal Jg, 
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but it does not tell all invariants of Iq- To be precise, although Theorem 12.11 shows that 



M = Lift(Fi) U Lift(F2) U Quad(G'i,G2) 

generates the cut ideal Jg, when Fi and F2 are minimal generating sets of Iq^ and /ca, 
the set M need not generate minimally. This happens in Example 12.51 Furthermore, 
we do not know how taking toric fiber products affects the Cohen-Macaulay type. For 
instance, the usual Segre product of two Gorenstein ideals need not be Gorenstein. 

3.2. Conjectures. We now present some conjectures inspired by our computations. Our 
main observation is that many of the coarse invariants of the cut ideals seem to be pre- 
served under taking minors of the underlying graph. Recall that a graph if is a minor 
of G if if can be obtained from G by deleting and contracting edges. By the Robertson- 
Seymour Theorem on graph minors [T7], we may hope to characterize the class of graphs 
whose cut ideals satisfy some algebraic property by a finite list of excluded minors. 

The pro typical example of such a conjecture concerns the maximal degree of a binomial 
appearing in a minimal generating set of the cut ideal Iq- This number is ^i{Ig)- 

Conjecture 3.1. The set of graphs G such that fi^Ic) < k is minor- closed for any k. 

As evidence for Conjecture 13.11 note that two operations related to taking graph minors 
amount to taking faces of the corresponding cut polytopes. 

Lemma 3.2. (1) IfH is an induced subgraph ofG then Cut'~'(if) is a faceofCut°{G). 
(2) If H is obtained from G by contracting an edge then Cut^{H) is a face of Cut^{G) . 

Proof. For part (2), intersect Cut'~'(G) with the hyperplane Xij = where ij is the con- 
tracted edge. For part (1), intersect Cut'^(G) with the hyperplanes Xij = for all edges 
ij in G not incident to H, together with one extra condition Xij = for each connected 
component oi G\H, where ij is an edge incident to said connected component and H. □ 

This implies that generating degrees can only go down when passing to an induced 
subgraph or when contracting an edge: 

Corollary 3.3. (1) If H is an induced subgraph of G then ^{Ih) < f^i^c)- 
(2) If H is obtained from G by contracting an edge then fJ^^In) < f^i^c)- 

Proof. For any two toric ideals, we always have the inequality fi{Ii3) < f^i^A) whenever B 
is a face of A. Thus, the desired inequalities are a direct consequence of Lemma 13.21 □ 

Therefore, to prove Conjecture 13. H it would suffice to show that generating degrees are 
nonincreasing upon the deletion of edges. Note that the face property does not hold when 
deleting an edge, as seen by comparing Examples 11.11 and II. 2[ 

Conjecture 3.4. Let H be obtained from G by deleting an edge. Then ^i{Ih) < f^i^c)- 

The smallest instance of Conjecture 13. H namely k = 2 concerns those graphs G whose 
cut ideal Iq is generated by quadrics. We propose the following simple characterization: 

Conjecture 3.5. The cut ideal Iq is generated by quadrics if and only if G is free of K4 
minors (i.e. if and only if G is series-parallel). 
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If a graph G has Kn as a minor, then that minor can be reahzed by a sequence of edge 
contractions only. By Corollary 13.31 (2), the cut ideal of every graph with a minor 
has a minimal generator of degree 4. Thus, to prove Conjecture 13.51 we must show that 
graphs without K/^ minors have quadratically generated cut ideals. Graphs free of 
minors are known as series-parallel graphs. Every series-parallel graph can be built from 
K2 by successive series and parallel extensions. The series extensions are just 0-sums. 
Hence, to prove Conjecture 13. 5[ it would suffice to show that ^{Ig) does not increase 
when performing a parallel extension. 

Another conjecture, along the same lines as Conjecture 13.51 concerns quartic generators. 

Conjecture 3.6. The cut ideal Iq is generated in degree < 4: if and only if G is free of 
minors. 

In algebraic statistics, minimal generators of toric ideals are called Markov bases 0, [HI 
[22] . Thus, what Conjectures 13.11 13.41 13.51 and 13.61 are about is the complexity of Markov 
bases for moves among the N-valued functions on the cuts of a graph G. As we shall see 
in Sections 4 and 5, the underlying toric models [151 §1-2] are important in statistics, and 
this endows our computations and conjectures in this section with an applied relevance. 

From the more theoretical perspective of commutative algebra, it appears that Conjec- 
ture 13.61 also captures the class of graphs having normal and Cohen-Macaulay cut ideals. 

Conjecture 3.7. The semigroup algebra K.[q]/lG is normal if and only ifK.[q]/lG is 
Cohen-Macaulay if and only if G is free of minors. 

That K[g]//x5 is not normal and not Cohen-Macaulay can be seen in Table 1. The gap 
between the codimension (5) and the projective dimension (15) is remarkably large in this 
case (we note that the associated semigroup Q and its saturation Qsat differ by only one 
point). The property of being normal is preserved when passing from a semigroup algebra 
to a facial subalgebra. Hence we can deduce from Lemma [3.21 that every graph with a 
minor has a non-normal cut ring 'K[q\/lG- Thus, to prove a large part of Conjecture 13.71 
it would be sufficient to prove that graphs G which are free of minors have normal 
semigroup algebras K[g]//G'. Here we are using Hochster's Theorem, which states that 
normal implies Cohen-Macaulay among semigroup algebras [T2] . 

One question that remains is to characterize those i^s-free graphs G whose cut ideal Iq 
is Gorenstein. Being Gorenstein seems to depend in a complicated way on the structure of 
the graph G. In general, the Gorenstein property is not preserved under taking toric fiber 
products and, in particular, is not preserved under taking clique sums of graphs. We do 
not have a firm conjecture on the structure of those graphs whose cut ideal is Gorenstein. 

4. From Cut Ideals to Binary Graph Models 

We now explain the correspondence between certain cut ideals and the toric ideals of 
binary graph models. These are statistical models for 2 x 2 x • • ■ x 2-contingency tables, 
whose algebraic properties were studied by Develin and SuUivant in [B]. Our main result 
in this section (Theorem 14.11) states that binary graph models on n nodes coincide with 
cut ideals of those graphs on n -|- 1 nodes where one node is connected to all others. 
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Let G be a graph with vertex set V = [n] = {1,2, . . . ,n} and edge set E, and suppose 
that G has no isolated vertices. We introduce a polynomial ring with 2" unknowns, 

K[p] = K[pi^i^...i^ \ ii,i2,...,in^ {0,1}], 

and a polynomial ring with 4 ■ lE'l unknowns, 

K[b] = K[bt^ \^^J e{0,l},eeE]. 

The binary graph model is defined by the following homomorphism of polynomial rings: 

{k,l}eE 

The kernel of ipc is a toric ideal which we denote by Jq- The binary graph model of G is 
the zero set of Jq in P^"~^. In statistics, this toric variety corresponds to the hierarchical 
model for 2 x 2 x ■ ■ ■ x 2 contingency tables where the 2 x 2-margins on the edges of G 
are fixed. The Markov basis for this model consists of the minimal generators of Jo- 

The suspension of the graph G = {V, E) is the new graph G whose vertex set equals 
[n+1] = I^U{?7, + 1} and whose edge set equals EU {{i,n+l} \i G V}. Given any binary 
string i = ziZ2 ■ ■ G {0, 1}", we define the associated partition 74(i)|i?(i) of [n + 1] by 
the condition k G B{i) if and only if i^. = 1. Similarly, if A\B is a partition of [n + 1], with 
n + 1 E A, we get a binary string i(y4|i?) by reversing this procedure. This specifies a 
natural bijection between the 2" unknowns pi in ]K[p] and the 2" unknowns qA\B in 

Theorem 4.1. Let 7 be the ring isomorphism K[p] ]K[g] defined by pi ^ (lA{i)\B{i)- 
Then, 

We note that this theorem is already known at the level of the underlying convex 
polytopes. This is the content of Chapter 5 of [7]. The polytope underlying the toric 
ideal Jq is the marginal polytope or covariance polytope of the graph G. It is isomorphic 
to the cut polytope of the suspension G under the covariance mapping, as explained in 
[3, §5.2]. The identification of Jq with Iq in Theorem 14. II lifts the covariance mapping to 
the setting of toric algebra. Before presenting the proof, we discuss a few examples. 

Example 4.2. Let G = be the complete graph on three nodes. The homomorphism 

i)G takes the polynomial ring IK[pooo,Pooi,Poio,Poii,Pioo,Pioi,Piio,Piii ] to the polyno- 
mial ring 612, 512^ 512^ 513^ 513, 513^ 513^ 523^^^^ by sending p,,, to b^bH^. 

The kernel Jq is the principal ideal generated by poooPoiiPioiPiio — PooiPoioPiooPiii- The 
isomorphism 7 sends pooo ^ gi234|, Pooi ^ «?i24|3, Poio ^ gi34|2, Poii ^ gi4|23, Pioo ^ 
gi|234, Pioi ^ '?i3|24, Piio ^ gi2|34, Pill ^ gi23|4- The image of Jk^ under 7 is the principal 
ideal I^^ which is discussed in Example 11.11 Note that is the suspension of K^. □ 

Example 4.3. Theorem 14. II explains some of the coincidences between rows in our Table 
1 and the table on page 447 of O §2]. For instance, the ideal Jka — Ir^, is minimally 
generated by 20 quartics and 40 sextics. Or, if G is the edge graph of the bipyramid, 
denoted BP in [6], then its suspension G is the graph G207 in our Table 1, and the ideal 
Jbp = Ig207 is minimally generated by eight quadrics, 436 quartics and 2872 sextics. □ 
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The results in Section 3 of [6] imply the following corollary for cut ideals. Note that it 
is consistent with Conjecture 13.61 because the relevant suspensions G have no minors. 

Corollary 4.4. Let G he a cycle Gn or a complete bipartite graph K2^n- Then the cut 
ideal Iq of the suspension G is generated by binomials having degrees 2 and 4. 

The results in Section 4 of [6] provide counterexamples to a conjecture that seems to be 
implied by Table 1; namely, there exist graphs whose cut ideals have minimal generators 
of odd degree. The smallest such example for a binary graph model concerns the graph 
G = K2X K3, the edge graph of the triangular prism, whose graph ideal Jq has a minimal 
generator of degree 3. The suspension of this graph, which has seven vertices, has a cut 
ideal with an odd degree minimal generator. 



Proof of Theorem 4^.1. It suffices to show that there are a pair of homomorphisms a 



-> K[s, t] and [3 : ]K[s, t] ]K[6] such that 0g o 7 = a o ipQ and ipc ° 7~"^ = j3 o 0g. 
The maps a and /3, restricted to ^[p\/Jg and K[g]/Jg respectively, will then lift to the 
isomorphism 7. To do this correctly, we extend K[s, t] and ]K[6] to allow fractional powers 
of the unknowns. Which fractional powers are needed will be clear from the context. 
We define the map a : K[6] ^ t] as follows: 

1 1 1 1 

Vkl J, J. J.deg(fe) ,dSi(!) ikl ,dcg{fc) dcg{0 



h'^^ i-^ a „dcg(fe) , dcg(0 ,kl , dcg(fe) d5i(Ty 

WlO bkli'k,n+l''l,n+l^ "U ^ ''kl^k,n+l^l,n+l- 

Here deg(/c) denotes the degree of the node k in the graph G, and similarly for the node /. 

We wish to show that a satisfies 0gO7 = aotpQ. To do this, we look at which unknowns 
Skh^ki appear to which powers in the monomials cx{'>Pg{P\)) and 0g(7(pi)). An unknown 
Ski appears in a(V'G(Pi)) with multiplicity one if and only if i^ii G {01, 10} if and only 
if {k,l} G Cut(y4(i)|i?(i)) if and only if Ski appears in 0g(7(pi)) with multiphcity one. 
A similar argument shows that tki appears with the same multiplicity in both a{il)G{P\)) 
and 0g(7(pi)). To check the multiplicity of Sk,n+i (and similarly for tk^n+i), note that 
the fractional powers guarantee that Sk,n+i appears in a{il)G{P\)) if and only if it has 
multiplicity one in oi{iIjg{p\))- This happens if and only if 2^ = 1 if and only if (/c, n + 1) G 
Cut(A(i)|5(i)) if and only if Ski appears in 0g(7(pi)) with multiplicity one. 

We now define our second ring homomorphism (3 : ]K[s,t] — ^ ]K[6] as follows: 

SM+i^ n Kh'^iY~'-B, tk,n+i^ n WiY~'-B 

l:{k,l}£E i.{k,l}GE 

ski^{b'o[b1l,)K hi ^ {b^l,b'A)K 

Here B denotes the product of all unknowns in K[6] raised to the power l/2n: 

n n C'S)'^- 

{k,l}£Ei,j£{0,l} 

To prove that f3 satisfies ipG 07"^ = /5o0gwe compare the multiplicity of b^J in 
i^G{l~^{QA\B)) and P{(pQ{qA\B))- By symmetry, it suffices to analyze the case ij = 00. 
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For fixed k,l, the unknown 6qq has muhiphcity one in ipG{'y~^{(lA\B)) if and only if 
{k,l} ^ Cut{A\B) and k,l E A. Here, 6q^, b^^, occur with multiphcity zero. 

Now we analyze the multiplicity of in P{4>Q{qA\B))- Suppose we are in the case 
{k,l} ^ C\it{A\B) and k,l E A. This means that tkitk,n+iti,n+i is a factor of 
Looking at the expansion of (3{4>Q{qA\B)), aside from the factor B^, the only multiplicands 
which possibly contain 6qq are tki, tfc,n+i, ti,n+i- The first contributes (&oo)^! second and 
third contribute nothing, and the factor of 5" contribute (600)2 for a grand total of ^qq- 
On the other hand, 601 appears with multiplicity zero because t^i and t^^n+i contribute 
nothing, t/^n+i contributes (^oi)"^' -B" contributes (690)^ • ^ similar argument shows 
that b^Q and also appear with multiplicity zero. This agrees with the multiplicity of 
b^j in iPg{i~^{(1a\b))- This completes the proof of Theorem 14. 1[ □ 

5. From Jukes-Cantor Phylogenetic Models to Cut Ideals 

In this section we apply cut ideals to phylogenetics. Our main result (Theorem 15.51) 
states that cut ideals of graphs with n nodes are precisely the binary Jukes-Cantor models 
on cyclic split systems on n taxa. This class includes the Jukes-Cantor models on phy- 
logenetic trees whose algebraic properties were studied in [3] and [19]. We rederive the 
quadratic Grobner basis for these ideals by relating Theorem 12.11 to [T9l Theorem 21]. 

The extension of statistical models of evolution from phylogenetic trees to split systems 
is due to David Bryant, who described these models in [2]. This extension has the double 
advantage of being useful for biological applications and leading to a richer mathematical 
theory. In what follows we give an algebraic introduction to Jukes-Cantor models for 
arbitrary split systems. Later on we specialize to split systems which are cyclic, and 
hence most relevant for the NeighborNet method [3]. This will take us back to cut ideals. 

5.1. The one-parameter model associated with a single split. We consider a set 
of n taxa labeled by [n] = {1, 2, . . . , n}. Each Jukes-Cantor model is a subvariety of the 
(2" — l)-dimensional projective space p2"-i whose coordinates we denote by Pi^-.-i^- The 
coordinate Pi^.-.i^ represents the probability of observing the states ii, . . . G {0, 1} at 
the taxa. We shall employ a linear change of coordinates which is known as the Fourier 
transform or Hadamard conjugation^ see [T^ §4.4] and [121 §2]. The Fourier coordinates 
are here denoted fji...j„, and they are related to the probability coordinates as follows: 



where the sum is over all elements {ii, . . . ,in) of the abelian group (Z/2Z)"'. It is very 
easy to invert this linear transformation. Namely, we have 



where the sum is over (ji, . . . , j„) G (Z/2Z)". 

A split {C, D} is a partition C U D = {1, . . . ,n} of the set of taxa such that n E D. 
We fix a split {C, D} and we introduce one free parameter u. In statistical applications, 
this parameter u would range over real numbers between and ^. In algebraic geometry 
we allow any point (mq : Mi) on the complex projective line P^, with Uq = 1 and Ui = u. 

We map the w-line into the probability space p2"-i by setting 



(5.1) 




(5.2) 
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• fji-ju = if ji H h j„ is odd. 

• fji-jn = uo if 'Ek&cJk and J^keaJk are both even. 

• fji-jn = ui if Y^kdcik and ^k&oJk are both odd. 

This hne in p2"-i is the Jukes-Cantor model associated with the spht {C, D}. Using the 
transformation flS.ip . we can express the parameterization in probabihty coordinates: 

• Ph-in = {uo + ui)/A if zi = ■■■ = «„. 

• Pii-in = (^0 ~ '"i)/4 if ^fc = 1 for all A; G C and = for all I G D. 

• Pii-in = (^0 ^ ^i)/4 if Zfc = for all A; G C and ii = 1 for all / G D. 

• Pij^...i„ = in all other cases. 

In summary, the Jukes-Cantor model for a single split is a straight line in P^"~^. Given 
two points in this model, we can multiply their Fourier coordinates, one coordinate at a 
time, and we get a new point in the model. Thus the model is a semigroup with respect 
to multiplication of Fourier coordinates. The model is a line which is also a toric curve. 

5.2. The Jukes-Cantor model defined by an arbitrary split system. A split system 
is simply a collection of r distinct splits of [n] = {1, . . . for some positive integer r: 

S = {{C,,D^},{C2,D2},...,{Cr,Dr}}. 

Each split {Cj, Di} specifies a one-parameter Jukes-Cantor model, which is a semigroup 
under multiplication of Fourier coordinates. We define the Jukes- Cantor model of S to 
be the semigroup generated by the r one-parameter models of the splits {Cj, Di} G S. 

Explicitly, the parametrization of this Jukes-Cantor model is given as follows. The pa- 
rameter space is the direct product of r copies of the projective line P^. The homogeneous 
coordinates of the i-th projective line are denoted [uq : u\). There are precisely 2"~^ 
nonzero Fourier coordinates fj^---j„- They are indexed by the group 

(Z/2Z)^,,, = { (ji, . . . , e (Z/2Z)" : + . . . + j„ is even }. 

Each nonzero Fourier coordinate is expressed as a monomial of degree r in the parameters: 

(5-3) fn..,. = n -L.c,.V 

Since this parametrization is given by monomials, the ideal of algebraic invariants of the 
Jukes-Cantor model is a toric ideal in the Fourier coordinates. This toric ideal is the 
kernel of the ring map fl5.3p and we denote it by JC^. It lives in the polynomial ring 
K[/] whose generators are the 2"^-*^ Fourier coordinates fji--j„ indexed by {Z/2Z)^^^^. 

It is important to understand that Jukes-Cantor models are toric varieties, since JCs is 
a toric ideal in the Fourier coordinates, but Jukes-Cantor models are not toric models (i.e. 
log-linear models or discrete exponential families) in the sense of [T5l §1.2] because JCs 
is not a toric ideal when rewritten in the probability coordinates Pi^...i„ via the Fourier 
transform (15.11) . 

Proposition 5.1. //S consists of r splits then the Jukes-Cantor model is r- dimensional. 

Proof. We can write the 2"~^ nonzero monomials in the parametrization (15.31) as the 
columns of a zero-one matrix A with 2r rows, one for each unknown Uq and u\, as in [18] 
or in [151 §l-2]- The rows of this matrix span an r -|- 1-dimensional linear space. This 
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implies that the semigroup algebra IK[/]/JCs has KruU dimension r + 1, and hence the 
associated projective variety (which is our Jukes-Cantor model) has dimension r. □ 

Jukes-Cantor models for split systems do indeed generalize the familiar models associ- 
ated with trees. Let T be a tree with leaves labeled by [n]. Every edge of T defines a split 
{C, D} of [n]. We write E(T) for the set of splits coming from all the edges of T. 

Remark 5.2. JCs(t) equals the usual Jukes-Cantor model associated with the tree T. 

Proof. This is seen by comparing the parametrization for split systems in (15 .Sp with that 
given in [19], §3] for group based models on trees. The condition that J^ked 3^ is even 
in the split system representation is replaced with the condition that XlfcGA(e) is even 
where A(e) is the set of leaves below the edge e. The concept of being a "leaf below an 
edge" is equivalent to being on one side of a split. □ 

5.3. Cyclic Split Systems. We now turn our attention to the family of cyclic split 
systems. These split systems are particularly useful for representing and analyzing metric 
spaces in biology, as they can be drawn in the plane using NeighborNet [3]. 

Formally, we define cyclic split systems as follows. We draw a convex n-gon in the 
plane and label the vertices by 1, . . . ,n in clockwise order. Every line in the plane that 
does not pass through any of the vertices defines a split {C, D}. The complete cyclic split 
system S*^"') is the collection of all splits of [n] = {1, . . . , n} which arise in this manner. 

Remark 5.3. The number of non-trivial cyclic splits in S*^"^ equals n{n — l)/2. 

A cyclic split system is any subset of S*^"^. In other words, a split system E is cyclic if, 
for each split {C, D} G E, the set C is an interval of integers C = [k, I] = {k, k + 1, . . . ,1}. 

Now we will show that every cyclic split ideal JC-£ is a cut ideal. We associate with 
each cyclic split system E a graph Gy. with vertex set [n] as follows. For each cyclic split 
{C,D} G E where n E D and C = [k,l], we introduce the edge {k — 1,1} in Gy.. (Here 
:= n). Thus G^ is a graph with one edge for each split in E. The representation of a 
cyclic split system E by its graph Gj: is very natural as the following proposition shows. 

Proposition 5.4. Let T be a planar tree with leaves labeled cyclically 1, . . . ,n and E(T) 
the associated cyclic split system. Then the graph Gy{t) consists of the edges in the 
subdivision of the convex n-gon which is dual to the tree T. 

Proof. The proof of this result is straightforward. The idea is illustrated in Figure 1. □ 

We now come to the main result in this section. We define a bijection between the set 
of all 2"~^ cuts of [n] and the set (Z/2Z)"ypj^ of binary strings that sum to zero. If A\B is 
any cut then the corresponding binary string ji ^2 ■ ■ ■ jn is defined as follows: 

jk = 1 if {/c — 1, A;} G C\xt{A\B) and jk = otherwise 

It is easy to see that ji + ■ ■ ■ + jn is even and that the cut y4|i? is uniquely encoded in the 
string jij2 ■ ■ ■ jn- This bijection defines an isomorphism of polynomial rings r : K.[q] 
K[/] by sending the unknown qA\B to fji - j„- 

Theorem 5.5. Let be a cyclic split system and Gy, the associated graph. Then the 
Jukes-Cantor model JCy equals the image of the cut ideal Iq^ under the isomorphism r. 
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Figure 1. A tree with seven leaves and corresponding subdivision of the 7-gon 

Proof. To see that the preceding bijection between cut coordinates and Fourier coordinates 
gives an isomorphism between the cut model for Gs and the Jukes-Cantor model JCs 
we must define an appropriate bijection between the parameters. This bijection between 
parameters is induced by the map which sends a cyclic split in E to an edge of the graph 
Gy,. Namely, we identify the parameter space associated to the split {Gi,Di}, where 
Gi = [k + 1, /], with the parameter space associated to the edge {k, 1} in via 

(5.4) (4 : u\) = [tki : Ski). 

Now, the unknown Ski appears in the squarefree monomial </'Ge(?^|b) if ^^^d only if 
{fc, /} G G\xt{A\B) if and only if j^+i + ■ ■ ■ + ji is odd if and only if Yliui^d '^^^ 
and only if the unknown u\ appears in the squarefree monomial on the right hand side of 
(15. 3p . Likewise, tki appears in 0Ge(9a|b) if and only if Uq appears in the right hand side 
of (15.31) . This shows, modulo the identification (15. 4p . that the image of the cut coordinate 
qA\B under the map (pc^ equals the image of the Fourier coordinate fj^-..j„ under the map 
(15. Sp . Therefore, both maps have the same kernel, and we conclude JCs = t{Ig^). □ 

Example 5.6. Let S = be the complete cyclic split system on four taxa, i.e., 

(5.5) S = {{12, 34}, {23, 14}, {1,234}, {2, 134}, {3, 124}, {123, 4}}. 

The associated graph Gs is the complete graph on {1,2,3,4}. With the ordering of the 
splits as in (15.50 . the map r and the Jukes-Cantor parametrization (15.30 are given by 
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Under the identification (15 ■4p . this coincides with the parametrization in Example ll.li 
The Jukes-Cantor ideal for the complete split system on four taxa equals 

Jt. = ( /oooo/oioi/ioio/iiii ^ /ooii/oiio/iooi/iioo )• 

The ordering of the factors in the above monomials coincides with the lexicographic 
ordering of the edges of K^. If we set Uq = u\ = = u\ = 1 in the parametrization, 
then we get the 4-cycle in Example \1.2\ which represents the Jukes-Cantor model for the 
star tree. This model is the same as the rooted claw tree Ki^^ in [131 Example 14]. □ 

5.4. Algebraic invariants for Jukes-Cantor models on cyclic split systems. The 

polynomials in the ideal Js are known as algebraic invariants in phylogenetics. When 
expressed in terms of the coordinates Pi^-.-i^ via (15.11) . these polynomials are the algebraic 
relationships which hold among the joint probabilities for all distributions in the model. 
Using Theorem 15. 5[ we can now translate our results and conjectures about cut ideals to 
the setting of Jukes-Cantor models. We begin by giving a new proof of a known result. 

Corollary 5.7. p!9| Theorem 2 (a)] Consider the Jukes-Cantor model for any trivalent 
tree T with taxa [n] . Then the ideal JCy.{t) has a Grobner basis consisting of quadrics. 

Proof. By Proposition 15.41 and Theorem 15. 5[ we have JCj^ij) = Igi where G is the edge 
graph of a triangulation of the n-gon. Such a planar graph can be decomposed into 
triangles using 2-sums. The result hence follows from Theorem 12. 1[ □ 

We now discuss the Jukes-Cantor ideals JCe for some other cyclic split systems. Each 
of the graphs G in Table 1 corresponds to such a split system. Namely, for each edge {/c, /} 
of G we introduce the cychc split {C, D} where G = {k + 1, k + 2, . . . ,1} and D = [n]\G. 

The complete graph Kn corresponds to the complete split system E'-"'-' . Table 1 reveals 
that the algebraic invariants for S*^^-' are generated in degree < 6 and the algebraic invari- 
ants for E*^^-* are generated in degree < 10. Conjectures 13.51 and Conjecture 13.61 translate 
into conjectures as to which Jukes-Cantor ideals JG^, are generated by quadrics and which 
are generated by quartics. Whether the generating degree h^JGy) for a cyclic split system 
can only decrease upon removal of a split is still unknown, in light of Conjecture 13.41 

Huson and Bryant [H] have shown that cyclic split systems, even if they do not arise 
from trees, always have useful representations by phylogenetic networks. However, this 
representation is generally not unique [HI Figure 5]. These split networks on n taxa 
are thus in many-to-one correspondence, via Theorem 15.51 to graphs with n vertices, and 
our results here shed light on the algebraic invariants of the associated statistical model 
[2]. One concrete application of this correspondence to phylogenetics will be the exact 
computation of maximum likelihood parameters for splits models as described in [T^ §6] . 

Example 5.8. Let n = 6 and consider the bipartite graph 3 where the bipartition sep- 
arates {1, 3, 5} from {2, 4, 6}. The corresponding split system S consists of the six trivial 
splits {{z}, [6]\{z}} and the three non-trivial splits {123,456}, {234, 156} and {345, 126}. 
This is the smallest split system whose split network is not unique. It is depicted in [T^ 
Figure 5]. Using our Table 1 in Section 3, we see that the corresponding Jukes-Cantor 
ideal JCs is minimally generated by 63 quadrics and 72 quartics. The semigroup algebra 
]K[/]/JCe is also normal and hence Cohen- Macaulay, by Hochster's Theorem [12]. □ 
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